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ABSTRACT 

We investigate the relationship between the quasi-thermal baryon-related 
photosphere in relativistic outflows, and the internal shocks arising outside them, which 
out to a limiting radius may be able to create enough pairs to extend the optically 
thick region. Variable gamma-ray light curves are likely to arise outside this limiting 
pair- forming shock radius, while X-ray excess bursts may arise from shocks occurring 
below it; a possible relation to X-ray flashes is discussed. This model leads to a simple 
physical interpretation of the observational gamma-ray variability-luminosity relation. 

Subject headings: gamma-rays: bursts - shock waves - X-rays: bursts 

1. Introduction 

Gamma-ray burst (GRB) light curves at 7-ray energies are often highly variable, and generally 
this is attributed to internal shocks occurring at some radius > 10 12 - 10 14 cm from the center of 
a relativistic outflow produced by a violent collapse or merger, beyond the "photospheric" radius 
at which the flow becomes optically thin to scattering by electrons associated with the baryons 
entrained. 

This photosphere is a source of soft thermalized radiation, which may be observationally 
detectable in some GRB spectra, and may also result in inverse Compton cooling of the 
non-thermal electrons accelerated in the shocks occurring outside it, thereby enhancing a hard 
GeV non-thermal component at the expense of the usual MeV synchrotron component. 

At small enough radii, however, the shocks can create enough pairs to re-establish a second 
photosphere caused by the pairs, and it is only beyond a limiting radius that the shocks remain 
optically thin to pairs. The most favorable region for shocks producing highly variable gamma-ray 
light curves is above this radius, while shocks occurring below it would lead to a second source 
of less variable radiation (Kobayashi, Ryde &: MacFadyen 2002; see also Ramirez-Ruiz & 
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Lloyd- Ronning 2002; and Spada, Panaitescu & Meszaros 2000), which is also X-ray rich. Bursts 
dominated by either the baryonic photosphere or by pair-producing shocks may be identified with 
X-ray excess bursts (Preece et al. 1996), while the latter resemble at least the harder examples of 
the proposed X-ray flash (Heise et al. 2001) sub-class of GRB. 

The existence of a limiting pair-forming shock radius provides also a scenario which combines 
recent work on the interpretation of afterglow light-curve breaks in terms of a jet opening angle 
(Frail et al. 2001; Panaitescu & Kumar 2001; Piran et al. 2002) or a universal jet shape (Rossi, 
Lazzati & Rees 2002; Zhang & Meszaros 2002; Salmonson & Galama 2002), and work on the 
gamma-ray variability- luminosity relationship (Fenimore & Ramirez- Ruiz 2001; Reichart et 
al. 2001). Identifying the pair shock radius as an approximate boundary above which shocks lead 
to more strongly variable gamma-ray light curves and below which shocks result in smoother 
X-ray rich light curves, a phenomenological jet model leads to a simple physical explanation of the 
quantitative form of the variability-luminosity relationship. 



2. Baryonic Photospheres 

Consider a relativistic wind outflow where the bulk Lorentz factor has a mean dimensionless 
entropy rj = L Q /Mc 2 = W 2 r]2 which varies (Ar/ ~ 77) on timescales t v ranging from a minimum 
dynamical timescale up to the maximum burst (wind) duration t w , 10~ 3 s < t v ^ t w . The flow 
starts from a minimum radius r Q = ct Vjjn i n = 10 7 r o j cm, and the Lorentz factor accelerates as 
r oc r up to a coasting (or saturation) radius r c ~ r Tf, beyond which it coasts as T = Tf. For a 
simple wind, neglecting finite shell effects, Tf = min[r/, 77*] where the value r?* is a critical value of 
the dimensionless entropy given by (Meszaros & Rees 2000) 

V* ~ = (LoaT/Airm^ro) 1 / 4 ~ lO 3 ^^) 1 / 4 . (1) 

Here £ Pj0 is analogous to the definition of the compactness parameter but using the proton instead 
of the electron mass. The coasting Tf values follow from the criterion that the proton drag 
time must be longer than the expansion time for protons to start to coast. Below the baryonic 
photosphere protons are naturally coupled to radiation, but in the optically thin region above the 
photosphere, if this occurs in the accelerating regime, the protons can still coupled to radiation 
and continue to accelerate out to a radius beyond the photosphere. The comoving density in the 
(continuous) wind regime is n! = {L /Airr 2 m p C'r]T), and using the above behavior of T below and 
above the coasting radius, as well as the definition of the Thompson optical depth in a continuous 
wind tt — n'aTir/T) we find that the baryonic photosphere where tt = 1 in the wind (w) regime, 
due to electrons associated with baryons, is 

r ph,w = I rji /3 7]~ 1/3 for r < r r]; ^) 
where for a wind r c = r minfry, r/*] is the coasting radius beyond which T = Tf = constant. 
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The accelerating and coasting behavior is followed on average also if the outflow consists 
of shells of duration t v = A Q /c > r Q , separated by intervals which could similarly be of order 
~ t v (or a superposition of several such frequencies), leading to an oscillatory modulation of the 
linear and coasting behavior. Aside from such modulation, in the optically thick regime the 
Lorentz factor can never exceed T < rj. However taking into account the finite shell structure, 
in the optically thin regime the coasting Tf will differ for some rj from the values min [r?,?7*] 
discussed in the wind problem. Taking for simplicity shells resulting from the minimum variation 
timescale A cty ^nin 

= r Q , the evolving comoving width of the shell is A' = [r, r Q rj, rj V] 
and the comoving volume of the shell is V' = 4-7rr 2 A / = [4-7rr 3 , A.irrjr r 2 , 47r?7~ 1 r 3 ] when r is 
[< r r], > r Q r], > r Q rj 2 ]. The comoving particle density n' = (Lt v /rjm p c 2 V') and the Thomson 
depth tt = tiotA! = 1 define a baryonic photosphere in the discrete shell (ds) regime 

ph,ds = rfrj^ 1 ! 2 for both r < r rj and r > r r\ . (3) 

For a wind made up of shells of approximate duration t v ejected at intervals of order t v , at high 
rj the shells move fast enough that a photon arising in one shell never crosses more than that one 
shell. At lower rj, however, a light ray can cross many shells before escaping, and the appropriate 
expression for the photosphere approximates that of the wind equation (|2|). The criterion for the 
latter to be valid is that r ^ A rj 2 = r a rj 2 , and the transition occurs at rj = rjt given by 

Vt = r,i /5 = 2.5xW 2 (L 52 r-}) 1 / 5 (4) 

Thus one has for the baryonic photosphere 

/~ 3 = vh' 3 for rj < r\ t \ , . 

-1/2 5/2 _!/ 2 f . l D J 
] i = rj t rj 1 for rj > r/t, 




where the first (wind regime) occurs only in the coasting regime, while the second (shell regime) 
applies partly in the accelerating and partly in the coasting regimes. These regimes differ from 
those in Meszaros Sz Rees (2000) by having a break at rjt = rjj instead of rj*, and by having 
a slope -1/2 above rj t instead of -1/3 (due to the shell regime, neglected in our previous paper). 
The photosphere is in the coasting wind regime for rj < rjt, in the coasting shell regime for 
f]t ^ V — r ?* / ' 3 ) an d in the accelerating shell regime for rj > rjt^ 3 (see Fig. 1). When the photosphere 
is above the coasting radius, the final Lorentz factor is just ~ rj. When the photosphere 
is below the coasting radius, the baryons continue to be dragged by the radiation above the 
photosphere until t' drag ~ (m p c 2 / 'ccjtu'^) exceeds t' exp ~ r/cT, out to a radius rd rag /r < rj 2 
where T ~ min^r/ -1 / 2 , rj*}. The final Lorentz factor is thus Tf = [rj, rj^rj^ 1 ! 2 , rj*] for values of 
[rj < rjt^ 3 , rjt^ 3 < rj < rj 2 , rj > rj 2 ]. This results in shock radii (Figure 1) which are divided into 
three regimes (instead of the two in Meszaros &; Rees 2000, where in the wind regime the values 
rjt and rfj 3 were collapsed into a single rj*). 

In units of the initial total luminosity L and initial temperature at r Q , 
G D = kT /m e c 2 ~ 21}^ r^ 2 (i.e. T ~ 1 L^ 2 i r o y 2 MeV), the lab- frame baryonic photospheric 
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luminosity L p h and dimensionless temperature @ p h behave as 

L q ( {r P h/r c )~ 2/3 = (r//v*) 8/3 = Vt 2/3 (ri/Vt) 8/3 for 77 < r) t ; 

-j- = -ft- = { (r p/l /r c )- 2 / 3 = v* 1/3 (v/v*) = Vt 2/3 iv/Vt) for 7? t < 77 < f]t /3 ; (6) 

1 f 4 /3 

LI tor TJ > 77* . 

Thus (L ph /L ) = (Oph/Qo) = r/7 8/15 = 2.5 x 1(T 2 (T p/l ~ 25/(1 + z) keV in the observer 
frame) for rj = rj t = rjt^ ~ 250; (L ph /L ) = (@ p h/@o) = f]* 1 ^ = f° r V = »7* = 1° 3 ; and 
(L p h/L ) = {Qph/Qo) = 1 at 77 = r/*/ 3 ~ 10 4 (where the photosphere occurs at the coasting radius). 

The internal shocks, which occur in the coasting regime at radii r/r Q = (A /r )r/ 2 > rj 2 , 
produce a shock photon luminosity 

L s h = Ce^iLo ~ 10~ 1 e eil / 3 ej il / 4 L o , (7) 

where the shock efficiency e s h,-i = 10~ 1 e ej i/3ej 1 /4 is a bolometric radiative efficiency when the 
cooling timescale is shorter than the dynamical time. Similarly the magnetic luminosity (if the 
turbulent field energy es = (l/3)e B1 / 3 is in equipartition with that of randomized protons and 
electrons) is L B = e B eiL ~ 10~ 1 e Bil / 3 ej il / 4 L o ~ L sh . Thus, L ph < L sh ~ L B for rj < rj t ~ 250; 
Lph < L sh ~ L B for rj t < 77 < 77* ~ 10 3 ; and L p/l > ~ L B for t/* < 77 < t/* /3 ~ 10 4 . This 
means that for 77 > 77* ~ 10 3 the baryonic photospheric component dominates the non-thermal 
internal shock component in a bolometric sense. This will lead to inverse-Compton cooling of 
the non-thermal electrons accelerated in the shocks, causing a weakening and softening of the 
nonthermal synchrotron spectrum of the shock, at the expense of a hard GeV) inverse Compton 
component, while most of the energy will be in a thermal X-ray component. 

The BATSE 7-ray luminosity is broad-band in nature, and can be written as 
Lj ~ (1/5)(1 — eic)L s h ^ (l/5)L s h where ejc is the IC efficiency, with a peak synchrotron 
frequency depending on the comoving magnetic field value B' . For low values of 77, shocks occur 
closer in, leading to higher B' and harder synchrotron peaks. For rj ~ 77^ the baryonic thermal 
X-ray photosphere may be responsible for the X-ray excess BATSE bursts (Preece et al. 1996). 
For lower 77 the photospheric thermal peak is even softer, while the shocks occur closer in and 
produce harder synchrotron peaks approaching the upper, less sensitive end of the BATSE band, 
which could lead to an apparent dominance of the soft X-ray thermal photospheric peak. 

For higher values rjt ^ rj 77* the thermal peak tends to blend with the synchrotron peak, 
resembling the canonical non-thermal GRB spectrum, while for 77 ^ 77* a hard (^ MeV) thermal 
component would be predicted to dominate. 



3. Shocks Above the Pair-Radius and Variable 7-ray lightcurves 



When shells of mass (m/2) with Lorentz factors Ti and T2 collide, the mechanical efficiency for 
conversion of kinetic energy mc 2 (Fx +T2) into internal energy is e« = (ri + T2 — 2^TiT2)/(Ti +T2), 
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where as before we parameterize e« = (1/4)6^1/4. If a total of 2N shells are ejected which collide, 
and the total isotropic equivalent kinetic energy of outflow is Ei SO , the corresponding internal 
energy produced in the merger of two shells is E int ~ €iN~ l E iso . Of that, a fraction e e is given to 
electrons, and for a high radiation efficiency in the MeV range and a high compactness parameter 
(i.e. high efficiency of pair formation) a fraction of order 1/3 of the radiated energy could be 
converted into pairs, and the energy in pairs in the merged shell is 

E± ~ e e ei(l/3N)E iso ~ W^e^/se^/^N^E^ erg. (8) 

Assuming that T in in the range [r m ,rjvf], with T m <Tm, the observed radiation comes mainly 
from collisions involving shells at the extremes of this range and is maximized for T m <C Tm- 
Such merged shells move with a center of mass Lorentz factor T c ~ y/T m TM- For shells of initial 
lab-frame widths A D ~ r G , for radii above the shock radius r s h = r a r] 2 (which is also the "expansion 
radius" above which the comoving width tx r and the comoving volume oc r 3 ) the energy radiated 
in the shocks can be enough to create pairs which make the shocked shells optically thick to 
Thomson scattering, if rj is below a certain value for which the comoving radiation compactness 
parameter £' ~ {Lot /m e c r) 2j 1 (Meszaros & Rees 2000). Earlier simulations involving randomly 
ejected shells and (baryonic) electron scattering in shocks have indicated a tendency for more 
variable light curves arising in more distant shocks (Panaitescu et al. 2000; Spada et al. 2000; 
Ramirez-Ruiz & Lloyd-Ronning 2002), as expected since closer in the scattering depth is larger. 
Similar results are obtained numerically when pair formation is included, e.g. Kobayashi et 
al. (2002). Here we pursue a simplified analytical description. For shocks at increasing radii, 
the pair comoving scattering depth of the shells eventually drops to unity, t'± ~ n'±aT(r /T c ) 
~ (E±ctt I '47rr 2 r ' c ) ~ 1 at a characteristic limiting pair-producing shock radius 

r± ~ {E±o T /Anm e c 2 T c ) 1 / 2 ~3x l^ie^e^/^ 1 E^) 1/2 {T M ,zT m , 2 )- l/i cm, (9) 

(Kobayashi et al. 2002) where T m = 10 2 r m) 2, Tm = 10 3 Fa,/,3- This is in the discrete shell regime, 
as opposed to the wind regime used Meszaros & Rees 2000; in general the shell regime pair density 
exceeds the wind regime pair density by a factor (t w /t v N)(T c /T m ) 2 > 1, as expected since the 
same kinetic energy density is concentrated in shells rather than smoothed out. At this radius 
both the comoving scattering time and the pair formation time as well as the comoving pair 
annihilation time (nj_o"Tc) _1 become equal to the comoving expansion time r±/cT c . 

Shells with Lorentz factors Tm and T m ejected from a starting radius r Q at time intervals 
t v = A G /c collide at a radius r s ^ ~ ct^T^ 2j r Q r/ 2 . If this shock radius is outside the limiting 
pair-shock radius r± given by equation pairs do not form in the shock, whereas in the opposite 
case an optically thick pair region does form in the merged shell, which expands until it reaches the 
radius r±. The shocks which occur outside the limiting pair-shock radius r± are those for which 
the corresponding shells started out from r Q with a minimum time difference t v > t v ±, where 

t v± ~ r ± /cT 2 c ~ 0.2(e eil/3 e iil/4 ^ 2 - 1 J E 54 ) 1/2 r- 5 / 4 r- 5 2 /4 s. (10) 
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If the shell ejection time differences t v have random realizations between the minimum and 
maximum values [t v>m , t«,Af] over the total duration of the burst outflow % <L t V) M, out of the N 
shells ejected there will be, on average, a fraction (1 — t v ±/t v ^) which will lead to shocks outside 
the pair-shock radius. For a high radiative efficiency, a fraction 0.5e 7; i/2 of which is taken to be in 
the gamma-ray range, the isotropic-equivalent gamma-ray fluence of the shocks above the limiting 
pair-shock radius is approximately 

E 1 ~ (l/2)e %1/2 e e €iE iso l-^ ~4x 10 52 e y ,i / 2€e,i/ 3 ei,i/4^54 f 1 - ) e rg. (11) 

\ £«,Af / \ l v,M J 

Here t v ±/t v>M ~2x 10~ 2 (e el /3e i)1 / 4 A r 2 " 1 ^54) 1/2 ^}r^ 4 r^ 4 < 1, with t v± < t V;M 4 where 
t;, = lOt;,^ s is the burst duration. In this simple model Ey represents the energy in the variable 
7-ray component of the burst, which arises above r± and has variability on timescales <; t v ±. For 
t v ± <C tj,, Ery is insensitive to T m , but for short bursts or for t v ± O.ltb there is a dependence of 

—5/4 

EL on t v ± oc T m . For small T m the typical pair-shock radius r± is further out, and the minimum 
variability timescale t v ± is longer, with a consequently smaller variable E^ (fewer shocks occur 
outside the more distant limiting pair-shock radius). Larger T m lead to smaller limiting pair-shock 
radii, shorter minimum variability timescales t v ± and larger isotropic equivalent £L. 



4. Pair-producing Shocks and X-ray Rich Component 

For shocks occurring at r s h < r±, i.e. below the limiting pair-shock radius where shocks can 
result in pair formation, the scattering optical depth of the shocked shells can become Tj_ <; 1 
(even when the shock is above the baryonic photosphere given by equations Bf| ). Pair formation 
causes the same amount of shock energy to be spread among a larger number of particles (new 
pairs) than in a purely baryonic outflow, and inverse Compton losses due to up-scattering of 
its own photons (Ghisellini & Celotti 1999) become important. For a pair of shells undergoing 
a shock at r < r±, it is expected that pair-production acts as a thermostat, and for comoving 
compactness parameters 10 ^ £' ^ 10 3 the comoving pair temperature is ~ 3 — 30 keV, with 
Tj_ ~ few (Svensson 1987). The scattering depth per shock due to pairs is unlikely to be much 
larger, because the scattering and the pair-formation cross sections are comparable, and unless 
dissipation and pair formation occurs uniformly throughout the entire volume, down-scattering of 
photons above the pair threshold rapidly leads to self-shielding (Ramirez- Ruiz et al. 2002). As a 
specific example, we take ~ 10 keV and Tj_ ~ 3 for one shock, producing a comoving spectrum 
peaked near hu' ~ 3kT±/ 1 t'\ ~ 3T!j_ io T '± 2 3 keV. Since at any time there may be more than one 
shock at r < r±, the photons might encounter more than one shell before escaping (e.g. Spada 
et al. , 2000), and would also undergo adiabatic cooling between the shells by a factor ~ r s h/r±. 
These two effects combined could lower the escaping photon energy by a factor roughly estimated 
as C ~ 0.2^o.2- For a CM bulk Lorentz factor T c = (rjy/r m ) 1//2 = 300r Cj 2.5, the observer-frame 
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pair-producing shock radiation peak is at 

hv x ,8h ~ 100 ri >10 T , ± 2 3Co. 2 r c , 2 .5[2/(l + z)] keV . (12) 

The peak energy ( |l2"l) is still substantially above the black-body value T s h,BB ~ 
4 (e e ,i/3 e i,i/4^54A r 2 _1 )" 1/8 rc,2{5 keV. The BATSE distribution of peak energies (Preece et 
al. 2000) has ~ 10% of bursts with hv pk < 100 keV, while the joint BATSE-BeppoSAX 
distribution of Kippen et al. (2002) shows that most X-ray flashes (XRF) have peak energies in 
the 20-100 keV range, with one exception at 3l 3 keV. A nominal value of hv pk ~ 20 keV can be 
obtained from equation fll2] ) with, e.g. r c [2/(l + z)\ ~ 60 



For completeness, we note that in the extreme case where pairs are produced uniformly 
throughout the entire volume, thermalization and an equilibrium pair optical depth t± oc i' 1 ! 2 
might be achieved, where £' is the comoving compactness (Guilbert, Fabian & Rees 1984; Svensson 
1987), although we expect t± in this case to be much smaller. 

The energy in the X-ray component from shocks arising below the limiting pair-shock radius 
r±, integrated over the burst duration tb, is the complement of the 7-ray energy produced in 



shocks arising above r± [c.f. equation (11)]. The X-ray isotropic equivalent fluence is 



E x ~ (l/2)e k e e eiE iso ^ ~4x 10 5 %, 1/2 (e e)1/3 e ijl/4 £54) 3/2 iV 2 ^{T^T^ erg, (13) 

\tv,M J 

where t v ±/t Vt M is given below equation (|i~T|), = (1/2)6^1/2 is an efficiency factor to account for 
a fraction of order unity of the luminosity below r± which is re-converted into kinetic energy. This 
X-ray component could account for most of the harder X-ray flashes (Heise et al. 2002; Kippen 
et al. 2002), with hv p k ^ 20 keV, but if more XRFs are observed with peaks as low as 3-5 keV 
this may require additional considerations. On the other hand, the X-ray excess GRB discussed 
by Preece et al. (1996) have characteristics which, as a class, are close to those of the r s h ^ r± 
pair-producing shocks discussed in this section. 

The radiation from pair-shocks with r± ~ few would be subject to a moderate amount of 
time-smoothing At var ~ At var ,origT±, which partially degrades the original variability implied by 
the random ejection and shocking of shells. The smoothing would be more appreciable at the 
shorter timescales, where it would lead to a filling in of the narrow troughs between peaks (see 
also Panaitescu et al. 2000; Kobayashi et al. 2002 and Ramirez-Ruiz & Lloyd-Ronning 2002). This 
smoothing, however, would not be expected to affect the coarser time structure of the light curve, 
since not many scatterings are incurred before the photons are advected with the flow. 



5. Variability Dependence on 7-ray Luminosity 



An observational correlation (Fenimore & Ramirez-Ruiz 2001; Reichart, et al. 2001) has been 
reported between the isotropic equivalent luminosity L 7 and a variability measure V of the 7-ray 
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time profiles, of the form 

L 7 oc V 9 , where g ~ 3.3±?;f . (14) 

The operational definition of V is related to the normalized variance, or the root mean square 
of the deviations from a smoothed light curve. Observations of afterglows with breaks in the 
light curves are believed to indicate the presence of a collimated jet-like outflow. The simplest 
interpretation assumes a uniform jet cross-section (independent of angle out to a jet edge 9), in 
which case the variety of break times indicates a variety of jet opening angles and the data indicate 
an isotropic equivalent fluence anti-correlation with jet opening angle 9 (Frail et al. 2001), of the 
form L 7i j so oc 6~ 2 . Alternatively, the same data can be interpreted in terms of a non-uniform 
(angle-dependent) cross section jet with a universal jet pattern given by the same functional 
relation between the energy output as a function of angle L^(9) oc 6~ 2 (Rossi et al. 2002; Zhang & 
Meszaros 2002; Salmonson & Galama 2002). In the latter case the data is interpreted as sampling 
different off-sets between the observer line of sight and the jet axis. Norris (2002) and Salmonson 
& Galama (2002) analyzing the time-lag effects (see below) in a larger sample and including 
redshift and luminosity function effects, argue for a somewhat steeper angular index of -5/2, so 

L 7 oc 9~ p , where p ~ 2 - 2.5 . (15) 

In the previous sections we used an isotropic outflow but our results continue to apply to the jet 
case as long as T exceeds the inverse of the jet opening angle 9. The model interprets the variable 
7-ray luminosity as that portion which arise from shocks above the limiting pair-shock radius r±, 
characterized by a minimum time variability given by equation (|!o|) , t v ± oc L 7 r m . This is 
based on equation (|ll|) relating i? 7i j so ~ E 1 to £54, and the assumption that the average mean 
duration and redshift differences are overshadowed by source- intrinsic variations in EL and T m , so 
that approximately L 7 oc E y , and assuming that t v ± <C t Vt M- The crucial dependence of t v ± is 
through T m , rather than Tm, since it is r" m which determines the shock radius. It is reasonable 
to make the ansatz T m oc 9~ g , and a value q ~ 2 (e.g. MacFadyen & Woosley 1999; Kobayashi 
et al. 2002) follows from momentum conservation in a "sharp boundary jet" model where the 
energy and T are constant throughout its cross section but there is a range of opening angles 
(e.g. Frail et al. 2001). In a "universal jet profile" model where L and T vary as function of 9 
(Rossi et al. 2002; Zhang & Meszaros 2002; Salmonson & Galama 2002), a value q ~ 2 is also 
expected, e.g. if the baryon loading in the jet are approximately independent of 9 but the energy 
varies as 9~ 2 (e.g equation |]i~5]| ). Setting T m oc 9~ q , we have then t v ± oc L^j 2 Ym^ oc L 7 2p 5g ^ 4p . 
The variability V of the gamma-ray light curves could be expected to scale, in an approximate 
way, inversely proportional to a power of the minimum variability timescale, V oc t~ k min . An 
approximate argument shows that such an anticorrelation exists in the GRB data, with an index 
k ~ 2/3 (e.g. Ioka & Nakamura 2001; Plaga 2001). Identifying t v ^ m i n with t v ±, we have 

L 7 OC y4 P /[fe(5g-2p)] _ (lg) 

If one takes p = q, which may be too idealized, the theoretical relation is L 7 oc V 2 , which is 
comparable to the lower limit fit of Fenimore & Ramirez-Ruiz (2001); the same result is obtained 
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for p = 5/2, q = 2, k = 1. Using the nominal values p = 5/2, q = 2, k = 2/3 we get L 7 oc V 3 , in 
good agreement with the observed best-fit relation L 7 oc V 3 3 of Fenimore & Ramirez- Ruiz (2001). 

6. Discussion 

We have discussed the properties of the quasi-thermal baryonic photospheric radiation 
component in GRB. At high isotropic equivalent luminosities, this component can dominate 
the non-thermal shock component, and appears in the hard X-ray range in the source frame. 
Such sources may be identified with the X-ray excess (Preece et al. 1996) class of bursts. This 
photospheric quasi-thermal component can inverse-Compton cool the non-thermal electrons in the 
shocks above it, suppressing the MeV synchrotron component and enhancing an inverse-Compton 
GeV non-thermal component. For high dimensionless entropy r] = L/Mc 2 and low (z ^ 1) 
redshifts the quasi-thermal component appears at hard X-rays (and in extreme cases at 7-rays), 
whereas for high redshifts (and/or low rj) it appears at soft X-rays. We also have identified 
a new regime in the description of baryonic photospheres from relativistic outflows, which is 
valid at moderate to high rj. The value of the final coasting Lorentz factor of the outflow is 
not automatically the value it has when the flow becomes optically thin, and has three different 
possible values 77, r/^r/ — 1 / 2 , 77*, as discussed below equation (|j|), depending on the value of the 
initial dimensionless entropy rj. 

We have quantified the location of the outermost radius at which pairs can form in internal 
shocks, and have argued that highly variable gamma-ray light curves arise mostly from shocks 
above this limiting pair-shock radius. The pair-shock radius determines the approximate ratio 
of the fluences in a variable gamma-ray non-thermal component and in a less variable softer 
(^ 20-25 keV) X-ray component. The latter could also be responsible for X-ray excess GRB, 
and, for moderately low bulk Lorentz factors or moderately high redshifts T[2/(l + z)\ <L 60, 
would be similar to most of the currently known X-ray flash (XRF) bursts (Heise et al. 2001; 
Kippen et al. 2001), but additional considerations may be needed to fit naturally the softest (3-5 
keV) XRFs. Smoother X-ray components are also obtained from closer-in shocks neglecting pair 
formation (e.g. Ramirez- Ruiz & Lloyd- Ronning 2002; Spada et al. 2000), but smoothing and 
softening is stronger when there is pair formation (see also Kobayashi et al. 2002). This pair X-ray 
component is generally softer than that of the baryonic photosphere. When present, the pair 
photosphere enshrouds the baryonic photosphere, but its modest opacity t± ^ 3 is not sufficient 
to alter significantly the spectrum of the baryonic photosphere. One or both of these X-ray rich 
components may be present, depending on the bulk Lorentz factor and isotropic equivalent total 
energy of the burst, and criteria are discussed for the non-thermal 7-ray components to dominate 
over, or be dominated by, these X-ray components. 

An individual burst is characterized in Figure 1 by an average r\ = (L/Mc 2 ). The shock 
radius r s /j j0 plotted in Figure 1 is for the minimum variability time t v = t Q ~ 0.3 ms, and a second 
shock radius r 8 hz is shown for t v = 10 3 t o ~ 0.3 s. For t v ~ t Q and rj ^ r/t the corresponding shocks 
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occur below both the baryonic r p h and pair shock r± photospheres, leading to X-ray rich bursts 
whose variability is partially suppressed. For t v <; 10 3 t o and rj <; i] t shocks occur at or above 
both r p h and r± leading to hard 7-rays with large variabilty at <; 0.3 s. An individual burst may 
have several variability timescales present, leading to both types of components simultaneously. 
Preponderance of one or the other leads to a short timescale but low amplitude variability X-ray 
rich bursts or XRF, or to a classical hard GRB with large amplitude variability mostly at <; 0.3 
s. Roughly speaking, X-ray flashes would be expected from the region 77 < rjt ~ 250, and classical 
GRB from 77 > rjt. A baryonic photosphere component should be present at the begining of bursts 
and X-ray flashes, and in the troughs between harder peaks due to shock radiation. However, 
the farther beyond the coasting radius the photosphere occurs, the weaker its energy fraction is 
relative to the shock and/or e^ component, because its energy drops as r~ 2 / 3 . Low values of 77 lead 
to further-out, weaker baryon photospheres, and at the same time to harder, relatively stronger 
shocks ocurring closer in to the photosphere. If long variability timescales are absent and 77 <C r\t 
a soft baryon photosphere may be the most prominent component, but its total energy would be 
very low. A strong baryonic photosphere dominated burst (with quasi-thermal <; MeV spectrum) 
is possible (for shock efficiencies ^ 0.1) for 77 ^ i]t^ 3 , and such a component may be detectable 
already for 77 ^ 77* ~ 10 3 . On the other hand, slower 77 ^ rjt ~ 250 outflows are likelier to make 
X-ray rich bursts through a pair-shock component. 

We argue also that the relationship between variable gamma-ray radiation and the limiting 
pair-shock radius leads, using the phenomenologically inferred dependence between isotropic 
luminosity and jet angle, to a simple analytical interpretation for the observed variability- 
luminosity relation L 7 oc V s (e.g. Fenimore and Ramirez-Ruiz 2001; see also Kobayashi et 
al. 2002). The positive correlation between variability and harder vF v peaks discussed by 
Lloyd-Ronning & Ramirez-Ruiz (2002) also finds a qualitatively similar interpretation in terms of 
a higher variability corresponding to closer-in shocks, which are more specifically in the present 
model shocks occurring just above the limiting pair-forming shock radius. 

Several physical explanations have been proposed for the presence of a cutoff above about 
1 Hz in the power density spectrum of GRB light curves (Beloborodov et al. 1998) in terms of 
baryonic electron scattering (e.g. Panaitescu et al. 2000; Spada et al. 2000; Ramirez-Ruiz &: 
Lloyd-Ronning 2002). Here we point out a different explanation for this, which is based on the 
existence of a minimum 7-ray variability timescale (equation |l(|) above which shock radiation is 
free from smoothing by opacity from pair- formation. This is seen in Figure 1, which shows that 
shocks associated with the variability timescales t v ?Z 1 s and 77 ^ 150 common among observed 
GRB occur above the pair photosphere. 

This research is supported by NASA NAG5-9192, NAG5-9153 and the Royal Society. We are 
grateful to S. Kobayashi for valuable discussions. 
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Fig. 1. — Schematic plot of the baryonic photospheric radius r p h\ the internal shock radii r s ^ j0 , 3 
for two different variability timescales t a , 10 3 t o where t a = 0.3 ms; the thin shell pair-producing 
photospheric radius r±; and the thin shell coasting radius r c starting at r Q = ct Q , as a function of 
i] = L/Mc 2 . Distinguishing between the continuous wind and discrete shell regimes leads to a new 
characteristic value rj t = rjt and a distinct regime between 77* and r/^ 3 (equations (|],f§, see text). 
Shocks can only occur above r s h , and pair-forming shocks only occur below r±. Also shown is 
the external shock radius for two external densities n ex t = 1, 10~ 3 cm~ 3 , giving an upper limit for 
internal shock radii. 



